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Abstract 

The study of radiative (j) decays into scalar mesons, with subsequent decay into 
TTTT'j, constitutes an important topic at the electron-positron collider DAPHNE. The 
interference of the respective amplitude with the one for 7r"'"7r~7 production, where 
the photon originates from initial state radiation, will allow for unambiguous tests 
of models for the (j) 7/0 '^'^) amplitude. The forward-backward asymmetry of 
charged pions, which is a clear signal of such an interference, amounts up to 30 % 
and is at the same time quite sensitive to the details of the various amplitudes. The 
results for several characteristic cases are presented and their implementation into 
the Monte Carlo generator PHOKHARA is described. 
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1 Introduction 



The method of the radiative return, i.e. the study of events with initial state 
radiation (ISR) at electron-positron colliders, allows to measure the hadronic 
cross section over a wide range of energies in a single experiment, without 
changing the beam energy. Originally suggested already long time ago [1] , the 
idea has been revived in [2], after high luminosity electron-positron colliders, 
specifically 0- and 5-meson factories, came into operation. Already now the 
method has been used to study the pion form factor at DAPHNE [3] and the 
production of J/iIj{—>- fi'^fi~) [4] and of 7i~^tt~7i^ at BaBar [5]. It relies on the 
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factorization into a part describing the production of the virtual photon from 
the initial state (including radiative corrections) and a part describing the 
production of hadrons from the virtual photon, which essentially corresponds 
to the cross section for e+e" — > hadrons. This important feature is affected 
by contributions from the emission of photons from the final state (FSR). As 
discussed in [2] , suitable kinematic cuts can be chosen, which strongly reduce 
these effects in e^e~ 7r"'"7r~7. For B-meson factories, configurations with 
hadrons of low invariant mass in one hemisphere and a hard photon in the 
opposite hemisphere are completely dominated by ISR, with a small admixture 
from the "two step" process e+e~ 77*(— 7r+7r~7). The situation is more 
involved at the 0-factory DAPHNE. For photons at small angles with respect 
to the beam direction and centrally produced pions, ISR is again dominant. 
This is the kinematic region studied in [3] for a measurement of the pion form 
factor. In general, however, FSR will be more important at these low energies. 
Unavoidably, the input for any FSR analysis and Monte Carlo simulation will 
be model dependent and a careful experimental test of the model assumptions 
is required. The most recent and detailed analysis of this problem [6] has been 
performed with the Monte Carlo event generator PHOKHARA [7], where 
FSR for pions is modeled by scalar QED (sQED) (a different ansatz, leading 
to similar conclusions is presented in [8]). This ansatz is certainly adequate for 
most of the kinematic regions and for cms energies away from the resonance. 
For ^/s ~ 1.02 GeV, on top of the 0-resonance, an additional complication 
arises: the presence of the narrow (/)-resonance leads to an enhancement of 71717 
final states through radiative 0-decays into (dominantly scalar) resonances. 
This process is important in itself and has, for many years, been the subject 
of theoretical investigations (see [9,10,11,12,13,14,15] and references therein), 
where the decay chain (f) (/o(980) + /o(600))7 Trn'-f is generally assumed 
to dominate. First investigations, where this channel has been linked to the 
radiative return, have been published in [13]. Two distinct features allow to 
discriminate the radiative (f) decay from 77777 production through the radiative 
return on one hand and from FSR from charged pions on the other hand: 
Radiative decays to 7/0 will lead to 7r"'"7r~7 and 7r°7r°7 states with a relative 
weight 2:1; this is in contrast to FSR from pions, which, similarly to ISR, 
leads by construction to 7r"'"7r~7 only. FSR is, furthermore, strongly enhanced 
for photons coUinear to the pions. A second important feature, which will be 
the main subject of this paper, is the charge asymmetry. It arises from the 
interference between ISR with 7r"'"7r~ in an odd charge conjugate state, and 
FSR and radiative (f) decays, with 7r+7r~ in an even charge conjugate state. 
With the ISR amplitude being already fairly well determined, this interference 
will allow to pin down size and even the relative phase of the — > 77777 
amplitude in a model independent way. As a third possibility the dependence 
of all this distributions on the cms energy and the energy dependent relative 
importance of and continuum amplitude could be used. However, we will 
not dwell further on this last, fairly obvious, aspect. 



2 



2 Description of the model 

The subsequent analysis will not attempt to study the most general ampli- 
tude for the radiative decay. Instead we will use characteristic models with 
different resonance admixture and variable phases and demonstrate that these 
will lead to marked differences in the charge asymmetry. 

For definiteness, we have adopted two models describing (j) — > 7r7r7 decays (tt 
stands for both neutral and charged modes), used also in [13]: the "no struc- 
ture model" of [9] and the K^K~ model of [10]. The models were extended to 
include both /o(980) and /o(600) (a) intermediate states in the decay chain 
(j) (/o(980) + /o(600))7 — >• 7r7r7, an extension required to fit experiments 
[16,17]. More sophisticated version of the K^K" model, incorporating also 
the complete /o(980) -|- /o(600) mixing matrix, can be found in [11,12]. Con- 
tributions from the 7rp(^ 717) decay chain are not taken into account. 
Indeed KLOE data [16] show, that they are negligible for the 7r°7r°7 final state. 
The branching ratios BR(p^ — > 7r^7) and BR(p° — > ■n^'-^) are comparable [18], 
hence the same apphes to the 77+77^7 final state. 

At present this choice of models is not based on definite experimental informa- 
tion, which is still quite poor, but on their simplicity. To discriminate between 
different models, more experimental work is needed. However, as shown in the 
next section, a lot can be learned from studies of charge and charge induced 
asymmetries. 

The amplitude for FSR based on sQED, denoted A^sQed, is identical to the 
one used in PHOKHARA and is given by 



The NLO correction to FSR and the complete description of ISR, including 
NLO corrections, is also identical to the one used in PHOKHARA and de- 
scribed in [7,6,19]. The pion form factor is taken from [20]. 

For the sake of definiteness, we write down explicitely the amplitude Al<^ for 




+ (?2 + A; - qi) 



g2-e*(7) 
q2 ■ k 




(1) 



(/.^(/o(980) + /o(600))7^7r7r7 



(2) 
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which, when combined with production, describes the resonant process 
e+e" ^ 0* ^ (/o(980) + /o(600))7 rnr^ with the amphtude 



Here P (e(0)) and A; (e(7)) are and photon four momenta (polarization vec- 
tors) respectively. P = pi+p2 = Q + k in the case of e'^{pi)e~{p2) annihilation 
amplitudes, qi (52) is the tt"*" (tt") four momentum and Q = qi + q2- de- 
notes mass and its width. Unless stated otherwise, all numerical values 
of physical parameters are taken from [18]. The g^fyy coupling, deduced from 
r(0 e+e"), is g^../ = 7.765 ■ 10^^ GcV^. The phase «^ in Me+e-(0) becomes 
relevant for the interference between M^+e- (0) and both A^sQED and ISR am- 
plitudes and thus can be measured for charged pions. It is kept constant in 
the 'no structure' model. In the case of the K^K' model it is taken from [12]: 
— h^J s — 4m^, with h — lb° /GeV (thus also constant for fixed s). 

The function f(j){Q^) reads 



in the 'no structure' model, and 

94>K+K- T ( 



K X^V^KJ 

9fQK+K-9foTT-iT ia„ 9cTK+K-9amr 



(5) 



in the K^K~ model. The meaning of the three-particle couplings 9mim2m3 is 
self-explanatory and, assuming isospin symmetry, these couplings are identical 
for charged and neutral pion modes. M/^ and F/g (M^- and F^) are the mass and 
the width of /o(980) (/o(600)). Again for simplicity, we use constant widths 
for both resonances. The phase has to be substantially different from zero 
to fit the experimental data [16,17] (see below). /(•, •) is the K-loop function 
given in [10]. 

From the amplitude Ai^ (eq. 2) we have calculated the differential branching 
ratio dBr{(f) — > n'^n^'y) /dM^^j,, as measured in [16,17] {^/Q^ = M^,^). The 
parameters obtained in the fits to these data are shown in Table 1. The values 
of masses and widths obtained in the fits do depend on the choice of the model 
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Table 1 

The values of the parameters obtained by fits to the experimental data [16,17] 









'no structure' model 


model 




Mf^ [MeV] 




983.9 ± 1.4 


1003 ± 8 




T/o [MeV] 




35.4 ±3.2 


108 ± 7 




[MeV] 




588 ± 145 


519 ± 23 




[MeV] 




1653 ± 921 


319 ± 65 








(101 ± 11)° 


(98 ± 8)° 




9(i>foi9foTT-K 




2.84 ±0.17 






94'cr'y9o'n7T 




4.0 ± 1.5 




94,K+ 


K-9foK+K-9foTnT 


[GeV2] 




59.4 ±5.1 


9(j)K 


^K-9aK+K-9<TmT 


[GeV2] 




7.9 ±2.0 



and at this point it is difficult to interpret them as physical masses and widths 
of the resonances. 
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Fig. 1. Fits to the measured [16,17] differential branching rations: K^K 
model - solid line {x^ = 29/(53 d.o.f.)), 'no structure' model - dashed line 
(x' = 38/(53 d.o.f.)). 



Both models agree well with the data, as shown in Fig.l, and the destructive 
interference between /o(600) and /o(980) amplitudes is needed in both cases 
to describe the low M^^, tail of the distribution. As expected, the 7r°7r°7 
data alone cannot differentiate between different models. However, as shown 
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in the next section, charge asymmetries have an enormous analyzing power, 
and allow to test details of the radiative decays to charged pions. 



3 The analyzing power of charge asymmetries 



The importance of charge asymmetries for tests of the FSR model was dis- 
cussed in details in the context of scalar QED in [2,6]. At the 0-factory 
DAPHNE, the direct radiative 0-decays, not included in [2,6], play, however 
an important role [13]. These a priori small effects are enhanced by the reso- 
nant behavior of and their accurate analysis is compulsory. The interest in 
these asymmetries is twofold. On one hand they play an important role in the 
FSR 'background' estimate for the measurement of the pion form factor via 
the radiative return method. On the other hand, as we will show below, they 
allow for powerful tests of models describing radiative 0-decays. The asymme- 
try analysis is thus a source of rich experimental information complementary 
to the cross section measurement. 
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Fig. 2. Forward-backward asymmetry for different radiative cj) decay models com- 
pared with the asymmetry calculated witliin sQED (no /o) 



Let us start with the forward-backward asymmetry defined for tt 
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Fig. 3. Angular and dependence of charge asymmetries for different radiative (/) 
decay models compared with the asymmetry calculated within sQED (no /o), for 
centrally produced photons and pions 



, .^2^ _ N{9^+ > 90°) - iV(g,+ < 90°) . 2\ 

-^^^^^ > - N{9^, > 90°) + Ni9^, < 90°^ ) ' 



In this analysis the z- axis is chosen along initial positron direction. The cuts on 
the pion angles 20° < 6t^± < 160° are required for pions to be observed in the 
detector and the photon angular range 45° < 6^ < 135° was chosen to enhance 
the FSR effects. The results are shown in Fig. 2. The asymmetry is very 
sensitive to the relative phase between sQED and direct decay amplitudes 
(a^) and predictions of this phase [12] can be definitely tested. The assumption 

of a dependent phase from the complex loop function I i:;^, ^ j can also 

be tested, since the interference pattern can be measured as a function of 
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Q"^. The effects depend also on the photon angular range and more detailed 
studies are possible. Sizable effects are predicted and, as one can qualitatively 
see from preliminary KLOE results [21] on asymmetries, some of the models 
we use may already be excluded by the data. 

The angular and the dependence of the charge asymmetry for different 
models are shown in Fig. 3. Again large effects are observed, with marked 
differences between different models. The combination of cross section mea- 
surements for both charged and neutral pion final states with asymmetry mea- 
surements in the charged mode will allow for the most comprehensive study 
of the (f) 7r7r7 decay. 
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Fig. 4. Differential cross section and for different radiative 4> decay models compared 
with sQED (no /o) cross section (left plot) and their relative difference (right plot). 



The effects of the direct radiative decay on the differential cross section 
are shown in Fig. 4 for centrally produced photons and pions. They are big 
and the question arise if they can affect the KLOE [3] analysis of the pion 
form factor. The answer can be found in Fig. 5. With the models used in this 
paper one predicts effects up to 1% in the e~^e~ 7r"*"7r~7 cross section in 
the ''KK-model" and the "no- structure" model, if = 7r/2 is adopted, but 
significantly smaller values for a different choice of the phase (an event 
selection close to the one used by KLOE has been adopted). All models are 
compatible with KLOE data on ^ 7r'^7r''7 decay and thus the simultaneous 
analysis of ^ 7r'^7r°7 and e~^e~ 7r+7r"7 data is highly desirable. It will 
also allow for a full control of the FSR contributions to the e~^e~ 7r"'"7r~7 
cross section, hopefully reducing the error from FSR in a{e~^e~ —* tt^tt^), as 
measured via the radiative return, to a negligible level. 
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Fig. 5. The relative size of the direct radiative (j) decay on the differential cross 
section for an event selection close to the one used by KLOE [3]. 



4 Summary and Conclusions 



The interplay between radiative corrections to pion pair production in electron- 
positron annihilation on one hand and irn'y production through radiative 0- 
decays on the other hand has been investigated. The charge asymmetry is 
shown to be a unique signal for the interference between the amplitude for tttt 
production through the radiative return and tttt production from the radiative 
0-decay. This asymmetry can be studied as a function of the invariant mass of 
the TTTT-system and is particularly sensitive to the model for (j) — >■ tttt'j. This has 
been demonstrated by implementing model amplitudes for e^e~ ttttj 
into the PHOKHARA event generator and by studying the resulting distribu- 
tions. Indeed different models, which lead to seemingly indistinguishable tttt 
mass distributions in 0-decay, give rise to strikingly different charge asymme- 
tries. We conclude with a discussion of the impact of radiative 0-decays on 
the measurement of the pion form factor and demonstrate, that their effect 
can be kept well under control, if suitable kinematic regions are selected. 
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